Searching for tetraquarks on the lattice 
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We address the question whether the lightest scalar mesons a and K are tetraquarks. We 
present a search for possible light tetraquark states with J PC = ++ and 1 = 0, 1/2, 3/2, 2 
in the dynamical and the quenched lattice simulations using tetraquark interpolators. In 
all the channels, we unavoidably find lowest scattering states n(k)n( — k) or K(k)n( — k) 
with back-to-back momentum k = 0, 2ir/L, ■ ■ ■ . However, we find an additional light state 
in the 7 = and 1 = 1/2 channels, which may be related to the observed resonances a and 
K with a strong tetraquark component. In the exotic repulsive channels 1 = 2 and I = 3/2, 
where no resonance is observed, we find no light state in addition to the scattering states. 

It is still not known whether the lightest observed nonet of scalar mesons a, n, eto(980) and 
/o(980) [5] are conventional qq states or exotic tetraquark qqqq states. Tetraquark interpretation 
was proposed by Jaffe back in 1977 [T] and it is supported by many phenomenological studies, 
for example [3] . The tetraquarks, composed of a scalar diquark and anti-diquark, form a flavor 
nonet and are expected to be light. The observed ordering m K < m- ao (980) favors tetraquark 
interpretation since the 1 = 1 state [ds] [us] with additional valence pair ss is naturally heavier 
than the 1 = 1/2 state [sei] [dit] . 

It is important to determine whether QCD predicts any scalar tetraquark states below 1 
GeV from a first principle lattice QCD calculation. Previous lattice simulations [U [3] have not 
given the final answer yet. The strongest claim for a as tetraquark was obtained using the 
sequential Bayes method to extract the spectrum [4] and needs confirmation using a different 
method. Our new results, given in this proceeding, are presented with more details in [6j[7]. 

We calculate the energy spectrum of scalar tetraquark states with 1 = 0, 2, 1/2, 3/2 
in dynamical and quenched lattice simulations. Our dynamical simulation (a ~ 0.15 fm, V = 
16 3 x32) uses dynamical Chirally Improved u/d quarks [5] and it is the first dynamical simulation 
intended to study tetraquarks. The quenched simulation (a~0.20 fm, U = 16 3 x28) uses overlap 
fermions, which have exact chiral symmetry even at finite a. 

The energies of the lowest three physical states are extracted from the correlation functions 
Cij(t) = (O|e>j(i)0j(O)|O)- =3 = T,n Z ? z T e ~ En 1 with tetraquark interpolators O - qqqq, 
where Zf = (0|O, |?i). In all the channels we use three different interpolators that are products of 
two color-singlet currents [6] . In addition, we use two types of diquark anti-diquark interpolators 
in 1 = 0, 1/2 channels [5]. 

When calculating the 1 = 0, 1/2 correlation matrix, we neglect the so-called single and 
double disconnected quark contractions [5], as in all previous tetraquark studies. The resulting 
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states have only a qqqq Fock component in this approximation, while they would contain also 
a qq component if single disconnected contractions were taken into account [5]. Since we are 
searching for "pure" tctraquark states in this pioneering study, our approximation is physically 
motivated. 

All physical states n with given J c = ++ and I propagate between the source and the 
sink in the correlation functions. Besides possible tetraquark states, there are unavoidable 
contributions from scattering states 7r(fc)7r(— k) for I = 0, 2 and scattering states n(k)K(—k) 
for / = 1/2, 3/2. Scattering states have discrete momenta k — ^-j on the lattice of size L and 
energy (m 2 + fc 2 ) 1 / 2 + (m 2 K + fc 2 ) 1 / 2 in the non-interacting approximation. Our main question 
is whether we find some light state in addition the scattering states in J = 0, 1/2 channels. 
If such a state is found, it could be related to the resonances a or k with a strong tetraquark 
component. 

The energies E n are extracted from the correlation functions Cij (t) via the eigenvalues 
A n (t) oc e _£ ™(* _t °) of the generalized eigenvalue problem C(t)u a (t) = \ n (t,to)C(to)u n at some 
reference time to [S]- 
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Figure 1: The resulting spectrum E n for 1 = 0, 2, 1/2, 3/2 in the dynamical (left) and the quenched 
(right) simulations. Note that there are two states (black and red) close to each other in I = and 
I = 1/2 cases. The lines at I = 0, 2 present the energies of non-interacting 7r(fc)7r( — k) with k = j^- 
and j = 0, 1, \p2. Similarly, lines at I — 1/2, 3/2 present energies of ir(k)K(—k). 
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The resulting spectrum E n for all four isospins is shown in Fig. [TJ The lines present the 
energies of the scattering states in the non-interacting approximation. Our dynamical and 
quenched results arc in qualitative agreement. 

In the repulsive channel / = 2, where no resonance is expected, we indeed find only the 
candidates for the scattering states 7r(0)7r(0) and tt(^-)tt(— with no additional light state. 
The first excited state is higher than expected due to the smallness of 3 x 3 basis. Similar 
conclusion applies for the repulsive 7 = 3/2 channel with ttK scattering states. 

In the attractive channel I = we find two (orthogonal) states close to the threshold 
2m„ and another state consistent with n(^-)ir(— ^), so we do find an additional light state. 
This leads to a possible interpretation that one of the two light states is the scattering state 
7r(0)7r(0) and the other one corresponds to a resonance with strong tctraquark component 
(see more general discussion in [10]). In the attractive 7=1/2 channel we similarly find the 
candidates for the lowest two Tt(k)K(—k) scattering states and a candidate for a k resonance 
with a large tetraquark component. These results have to be confirmed by another independent 
lattice simulation before making firm conclusions. 

We investigate two criterion for distinguishing the one-particle (tetraquark) and two-particle 
(scattering) states in [7] . The first criteria is related to the time dependence of CV,- (i) and A™ (t) 
at finite temporal extent of the lattice. The second is related to the volume dependence of the 
couplings (0|Oj|n). 

The ultimate method to study a and k on the lattice would involve the study of the spec- 
trum and couplings in presence of qqqq O qq O vac <-> glue mixing, using interpolators that 
cover these Fock components. Such a study has to be done as a function of lattice size L in 
order to extract the resonance mass and width using the Luscher's finite volume method p~0l ITT] . 
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